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SUMMARY 


Values for stability and control parameters have been determined by use of 
the equation error method and the maximum likelihood method from maneuvering 
flight data for a low-wing, single-engine, general aviation airplane. The air- 
plane responses were excited from steady flights at different airspeeds using 
the stabilator, aileron, and rudder deflections. The model of the airplane is 
based on the equations of motion with the linear aerodynamics. From the repeated 
measurements, the two standard-deviation confidence intervals for the estimated 
parameters were established. These bounds are used for the comparison of param- 
eters determined by both methods and also for the assessment of an effect of dif- 
ferent input forms and power settings. The static parameters are also compared 
with results from steady flights. Using these comparisons, the best values of 
estimated parameters were determined and their accuracies specified. 


INTRODUCTION 

The National Aeronautics and Space Administration is currently involved in 
extensive general aviation stall-spin studies. During the research program, 
several airplanes have been tested in the wind tunnel and in flight, and more 
tests with other airplanes are anticipated. In undertaking the stall-spin 
research, the airplane dynamics in prestall regimes must be understood. For 
that reason part of the overall program includes the measurement of airplane 
transient maneuvers for the extraction of a complete set of stability and con- 
trol parameters. These parameters include aerodynamic derivatives and the val- 
ues of aerodynamic coefficients corresponding to steady flight conditions. 

There have been several previous attempts using systems identification to 
determine parameters of general aviation airplanes from unsteady measurements. 
These attempts differ in the amount of data available, estimation techniques, 
and verification of results obtained. In reference 1 the equation error method 
(regression analysis) is applied to measured longitudinal data corresponding 
to good excitation of the long- and short-period modes. The same technique 
is used in reference 2 for the determination of the lateral derivatives from 
flights with different values of thrust coefficients. The equation error method, 
based on a least-squares technique, is very attractive because of its simplicity. 
It can be easily applied to each of the equations of motion separately and pro- 
vides direct estimates of the unknown parameters. The resulting estimates are, 
however, biased as a consequence of the measurement errors in the input and out- 
put variables. 

A second procedure used in airplane parameter estimation is the output error 
method. Because it usually uses the maximum likelihood estimation, it is often 
called the maximum likelihood method. The airplane longitudinal and lateral aero- 
dynamic parameters obtained by this method are presented in references 3 and 4 
and are compared with aerodynamic derivatives obtained from wind-tunnel tests 
and theoretical predictions. The maximum likelihood estimates are theoretically 



superior to those obtained from the equation error method. These estimates are 
asymptotically unbiased, consistent, and efficient, provided that the model of 
an airplane is correct and the input variables are measured without errors. 
However, the maximum likelihood method applied to the problem mentioned is time 
consuming because of its iterative nature and because all equations of motion 
considered enter the estimation algorithm. In sane experiments, small variances 
of the measurement noise, unknown modeling errors, and a limited number of data 
points could substantially reduce the superiority of the maximum likelihood 
method to the equation error method. Under these conditions, both methods might 
provide identical values for the estimated parameters. Detailed description and 
comparison of both methods can be found in references 5 and 6. 

The purpose of this report is to document estimates of the stability and 
control parameters for one of the general aviation airplanes involved in the 
stall-spin program. The parameters are extracted from longitudinal and lateral 
maneuvers initiated from steady flights at different airspeeds. The airspeed 
range extends from the minimum airspeed at which the airplane can still be maneu- 
vered to the maximum airspeed in horizontal flight. The two methods already men- 
tioned were applied to measured flight data in an attempt to obtain more accurate 
values of the stability and control parameters for the test airplane. 

This report first describes the test airplane, instrumentation, flight tests, 
and data reduction. Then the mathematical model of the airplane is introduced, 
and the estimation methods are outlined. The results from both methods are then 
compared. The static parameters are also compared with the results obtained 
from steady flights. Last, the effect of input form and power setting in the 
estimated parameter values is demonstrated, the best values of parameters are 
determined, and their accuracies are specified. 


SYMBOLS 

A wing aspect ratio 

a l = 9c L,t/ 9a t 

ax/ayr^z reading of longitudinal, lateral, and vertical accelerometer, respec- 
tively, g units 

b wing span, m 

by constant bias error in variable y 

C D drag coefficient, D/qS 

Cl lift coefficient, L/qS 

C Lf t lift coefficient of tail, L t /qS 

Ci rolling- moment coefficient, M x /qSb 

pitching-moment coefficient, My/qSc 
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yawing-moment coefficient, M z /qSb 
thrust coefficient, T/qS 
longitudinal-force coefficient, Fx/qS 
lateral-force coefficient, Fy/qS 
vertical-force coefficient, Fz/qS 
wing mean aerodynamic chord, m 
drag, N 
= c h n o/^n 

forces along X, Y, and Z body axes, respectively, N 

terms in equations of motion defined by equations (AT 9) and (A20) 

function which represents state-equation model 

acceleration due to gravity, m/sec 2 

function which represents output-equation model 

sensitivity matrix 

stick-fixed center-of-gravity margin 

distance of center of gravity aft of leading edge of wing mean chord 
expressed in percent of c 

distance of aerodynamic center aft of leading edge of wing mean chord 
expressed in percent of c 

moment of inertia about X, Y, and Z body axes, respectively, kg-m 2 
product of inertia, kg-m 2 
cost function 

- nPT 


K<j, 


term defined by equation (B16) 


k<Sa'k<5r '^ZnB'^nZft terms defined by equations (B17) to (B19) 


L lift, N 

l distance of aerodynamic center of tail aft of aerodynamic center of 

airplane without tail, m 
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M 


m x ,m y ,m z 

m 


N 

-> 

n 


P 

p(z/0) 

q 

q e 


distance of aerodynamic center of tail aft of neutral point of 
airplane, m 

distance of aerodynamic center of tail aft center of gravity, m 
Fisher information matrix 

rolling, pitching, and yawing moments, respectively, N-m 
mass, kg 

main diagonal element of the M matrix 

number of data points 

measurement noise vector 

roll rate, rad/sec or deg/sec 

likelihood function 

pitch rate, rad/sec or deg/sec 

number of unknown parameters 


q 


R 


S 

S t 

s(y) 

s jj 

T 


t 

-> 

u 


dynamic pressure, -PV 2 , N/m 2 

2 

measurement noise covariance matrix 
yaw rate, rad/sec or deg/sec 
wing area, m 2 
tail area, m 2 

standard error of variable y 

main diagonal element of the M -1 matrix 

thrust, N 

time, sec 

input vector 


u,v,w 


V 


4 


velocity along X, Y, and Z body axes, respectively, m/sec 
airplane total velocity, m/sec 



N + •< + 


indicated airspeed, knots 

Vrp modified tail volume defined by equation (B3) 

x state vector 

Xa/Ya x- and y-coordinates of angle-of-attack vane relative to airplane 
center of gravity, m 

X(3,z3 x- and z-coordinates of sideslip vane relative to airplane center of 
gravity, m 

output vector 

measurement vector 

a angle of attack, rad or deg 

a v angle of attack measured by wind vane, rad or deg 

3 sideslip angle, rad or deg 

B v sideslip angle measured by wind vane, rad or deg 

S a aileron deflection (one-half of sum of left aileron deflection and 

right aileron deflection) , rad 

6 e stabilator deflection, rad or deg 

S r rudder deflection, rad or deg 

6 t trim tab deflection, rad or deg 

e downwash angle at tail, rad or deg 

0 unknown parameter 

0p predicted value of unknown parameter 

5 vector of unknown parameters 

0 pitch angle, rad or deg 

Xy scale factor error of variable y 

Vy residual of variable y 

p air density, kg/m^ 

a standard deviation 



4> bank angle, rad or deg 

4> yu phase-angle characteristics relating y and u variables, deg 

to angular frequency, rad/sec 


Aerodynamic derivatives {referenced to a system of body axes with the origin 
at the aircraft center of gravity, which is located at 20.6 percent of c) : 
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c m,o' Cnicx' Cmq' and c m<$ e defined in appendix A (egs. (A8) to (All)). 
Subscripts: 

E measured 

o trimmed condition 

t tail 

Superscripts: 

T transpose matrix 

-1 inverse matrix 

estimated values 
mean 

derivative with respect to time 

-V 

vector 

Abbreviations : 
c.g. center of gravity 

EE equation error 

ML maximum likelihood 

rms root mean square 


TEST AIRCRAFT AND INSTRUMENTATION SYSTEM 

For this study, a four-place, low-wing, single-engine airplane was used. 
The control surfaces included conventional ailerons, rudder, and all-movable 
tail (stabilator) . The basic geometric, mass, and inertia characteristics are 
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summarized in table I. The moments of inertia were measured for the airplane 
in its early test configuration. The airplane was later modified by the instal- 
lation of an onboard rocket system which is used primarily for spin recovery. 

The resulting changes in the airplane configuration affected only its mass and 
inertia characteristics. New moments of inertia were calculated from those pre- 
viously measured. 

An analog measurement system was installed in the airplane for recording 
control surface deflections, stick and rudder forces, airplane response vari- 
ables, and other quantities defining flight and engine conditions. Control posi- 
tion motions (input variables) were measured by rotary potentiometers directly 
attached to the control surfaces. An orthogonal triad of linear accelerometers 
was rigidly mounted on the center line of the cockpit floor at a location close 
to the allowable center-of-gravity range of the airplane (fig. 1). The sensitive 
axes of all accelerometers were alined to the reference axes of the airplane. 

Incidence angles were measured by a swiveling vane mounted on booms ahead 
of each wing tip (fig. 1). Because the corrected readings of both vanes gave 
identical results, only the angle-of-attack and angle-of-sideslip data from the 
right vane were used for the analysis. The indicated airspeed was obtained from 
the airplane's air data system which consisted of a simple total pressure orifice 
located on each side of the fuselage. Total temperature was measured by a sensor 
located on the top of the fuselage. The remainder of the instrumentation system 
included three rate gyros, attitude gyros, signal conditioning, power supplies, 
and tape recorder. These components were mounted on a rack behind the front 
seats as shown in figure 1. A summary of measured quantities used in this study, 
transducers, and static characteristics of corresponding channels is presented 
in table II. The root-mean- square (rms) errors were estimated from recorded 
signals during the preflight and postflight ground operation of the instrumenta- 
tion system with the airplane engine running. Both the resolution and the rms 
errors are referred to the digitized data. 

Table III presents dynamic characteristics of transducers used for the 
measurement of airplane response. These characteristics were obtained from 
dynamic calibration. The equivalent time constants given in the last column 
of table III represent the approximation of the transducer dynamics by a first- 
order system. 


FLIGHT TEST AND DATA REDUCTION 

Airplane responses were measured in six flights. Table IV summarizes per- 
tinent flight test conditions and the average mass and inertia characteristics 
of the airplane in these flights. Mass and inertia characteristics for each 
run analyzed were determined from the airplane take-off weight and estimated 
fuel consumption during the flight. 

The longitudinal and lateral modes were excited separately, primarily from 
the trimmed level flights at the airspeeds listed in table IV. For the investi- 
gation of power effect, perturbations were initiated from a steady climb with 
full power and from a steady descent with idle power. 
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In longitudinal flights, the inputs used were stabilator deflections having 
the form of a pulse, a doublet, or a combination of both. In the lateral case, 
both the rudder and aileron were applied simultaneously. Various forms of these 
inputs are shown later. In all cases, the a- and 3-traces were examined to 
determine that the atmospheric turbulence was negligible. 

The measured flight data were filtered with a 6-Hz low-pass filter and sam- 
pled at the rate of 20 samples per second. The sampled data were used to pro- 
duce automatic data tabulations, time history plots, and final tape for airplane 
parameter estimation. This tape included the following variables: time, true 

airspeed, incidence angles (right vane) , angular velocities, attitude angles, 
linear accelerations, control surface deflections, and incidence angles (left 
vane) . 

True airspeed was obtained from the indicated airspeed by applying correc- 
tions for measured position error of the static pressure system and by using 
the air density values computed from the measured air temperature and static 
pressure. The angle-of-attack vane readings were corrected for air upwash by 
a multiplication constant. This constant was estimated from steady horizontal 
flights by comparing longitudinal accelerometer and wind vane readings. The 
recorded linear accelerations were converted into the acceleration of the air- 
plane's center of gravity. The effective aileron deflection was computed as 
a mean value of the sum of the right and left aileron deflections. 

The next step preliminary to airplane parameter estimation included a 
compatibility check of measured response variables in steady and maneuvering 
flights. The relationship between variables a v , 0, az , and ax, and 4>, 

8 V , and ay was examined from the initial steady parts of various test runs. 
These data showed very small scatter in values of the longitudinal and lateral 
accelerations and sideslip. It was, therefore, assumed that the measurements 
of a x, ay, and 8 V were corrupted only by zero-mean random noise. Then the 
bias errors in ot v , 0, az» and 4> in the form of constant offsets were deter- 
mined. Similar bias errors in p, q, and r were found by assuming steady 
flight conditions. All these estimates were verified by the analysis of tran- 
sient maneuvers. The compatibility check of aircraft response variables in 
maneuvering flights included the prediction of V, B v , a v , 4> , and the esti- 
mation of constant bias errors in measured data. The technique used is based 
on airplane kinematic equations and an extended Kalman filter and is described 
in reference 7. 

Typical results from the compatibility checks are given in figures 2 to 5 
and in tables V and VI. In figure 2 the measured and predicted responses in 
V, a v , and 0 are compared. A similar comparison for the variables V, 8 V , 
a v , <t>, and 0 taken from one of the lateral maneuvers is presented in 
figure 4. The resulting residuals and the standard errors of the measured 
responses estimated from these residuals are included in figures 3 and 5. All 
the residuals, in general, indicate good agreement between measured and pre- 
dicted data. The larger discrepancy in measured and predicted airspeed is 
still within the rms error of the measuring instrument, which was- estimated 
to be equal to 0.89 m/sec. (See table II.) 
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Table V compares various estimates of bias errors in the longitudinal out- 
put variables. The error estimates from transient data were unaffected by the 
assumption that b ax = 0 and were close to those from the steady data. Because 

this pattern was observed in other maneuvers that were analyzed, only the esti- 
mates from steady-state data were used for corrections of the longitudinal 
responses. 

For the lateral case presented in table VI, the estimated bias error in 
ay is significantly different from the initially assumed zero. The inclusion 
of b a ^ as an unknown parameter in the model affected the estimate of b^ only. 

This particular run and other similar runs analyzed indicated some differences 
between the bias errors estimated from steady-state and transient data. These 
differences had, however, no significant effect on the estimated stability and 
control parameters. The errors estimated from steady-state data were, therefore, 
used for corrections of measured time histories of p, r, and 4>. 

MATHEMATICAL MODEL OF AIRPLANE 

A mathematical model of the airplane was formulated in the form of the state 
and output equations as 


= f(x,u,S,t) 


( 1 ) 


and 


■* 

y = g (x,u,0,t) 


( 2 ) 


where x, y, and u are the state, output, and input vectors, respectively, and 
® is the vector of unknown parameters. For the longitudinal motion, the four 
vectors in equations (1) and (2) have the form 


x T = [u,w,q,0] 

y T = [v,a v ,q,0 ,a x ,a z ] 

u T = 6e 

0T = jCx^'CXa^Z^'CZot^Zq'Czge^^'^'Cmq'CmSe] 


and for the lateral motion 
10 
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x T = [v,p,r,<j>] 
y T = [6 v ,p,r ,4>,a y ] 

uT = Pa'^r] 

® T = [ c Y,o' C '!£Q' C '!£ p ' c Y$ r ' C l ,o' C l$' C l p ' C l r ' C l$ a ' C ls r ' 
c n ,o' C ng /Cnp/Cn r ' c n 5a 'Cn<$ r ] 

The equations of motion are presented in detail in appendix A. Their form 
resulted from an examination of measured responses of the test airplane, and 
from wind-tunnel and flight-test results on similar airplanes. The coupling 
between the lateral and longitudinal motion during the measurement of lateral 
responses is included in the lateral equations by replacing the variables u, 
w, q, and 0 with their measured values. 


ESTIMATION METHODS 

The equation error (EE) method represents an application of regression 
analysis to each state equation separately; it is a method which minimizes the 
sum of squared errors satisfying the equation. The cost function for the state 
equation has the form 


N 

J r = ~ ^ [*rEi 
i=l 



where N is the number of data points and E denotes the measured quantity. 

The least squares solution for (§ r is obtained by finding the minimum of J r . 
The standard errors of the parameters are obtained from the information matrix 

in the so-called normal equations and from the residuals x rE i - f r i( x E ,u,0 r J . 
(See ref. 2.) ' 

In the airplane force equations the variables x rE were replaced by 
measured linear accelerations as indicated by equations (AT 2), (AT 3), and 
(A23) to (A25) . In the moment equations the angular accelerations were cal- 
culated from measured angular rates using splines. 

Detailed descriptions of the maximum likelihood (ML) method and of the 
computing program can be found in references 8 and 9. The ML method is based 
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on the maximalization of the log-likelihood function log p(z/3) . The comput- 
ing starts with the approximate values for the unknown parameters in the state 
equations and then iterates until the minimum of the cost function 


N 

i=l 


is found. In equation (4), z is the measured vector z = y + n, where n 
is the measurement noise vector. The measurement noise covariance matrix R 
is estimated as 



R 


1 

= diag - 
N 


N 



i=1 


(5) 


- 4 

where = z^ - yi are the residuals. 

After the kth iteration the new estimates of unknown parameters are found 
as 


®k = ©k-1 + A0 k 


( 6 ) 


$- 

where the vector A©^ is computed from the expression 


A ®k = — M ^ (0^_i ) 


) T 


(7) 


The matrix M is the Fisher information matrix. The inverse of the infor- 
mation matrix gives a lower bound on the error covariance matrix for the esti- 
mated parameters. Using the modified Newton-Raphson method, described in refer- 
ence 6, this matrix is approximated as 



In this expression H is the sensitivity matrix whose elements are 3y/9©j , 
where j = 1 , 2, . . . r qg and where q 0 is the number of unknown parameters. 
The main diagonal term of the information matrix mjj defines the sensitivity 
of measured output variables with respect to the parameter © j . (See ref. 10.) 

The main diagonal term of M -1 , s jj» is, therefore, the Cramer-Rao lower 

A 

bound on the variance of the parameter 0 j . It can be shown that 


s s 


m 


33 


33 


where the equal sign holds for the uncorrelated parameters. 


RESULTS AND DISCUSSION 
Longitudinal Characteristics 

The first part of the longitudinal data consists of 16 runs measured at 
different airspeeds. In these runs the transient motion of the airplane was 
excited from level flights by using the stabilator deflection in the form of 
a single pulse or doublet. Estimated parameters and their standard errors from 
one run at C L = 0.61 are summarized in table VII. All parameters included 
in the model indicate satisfactory identif iability as expressed by their stan- 
dard errors or lower bounds on these errors. Significant differences in the 
estimates using the EE method and ML method exist in the derivatives C Za 

and Cj^. From calculated parameter covariance matrices a strong correlation 

(the correlation coefficient about 0.94) was observed between the parameters 
C Z g and c z§ e and between C^ and Cmgg* regardless of the estimation 

method. 

Measured time histories and those computed by using parameters obtained 
by the EE method are given in figure 6 . Considerable disagreement is apparent 
only between measured and computed outputs V and ax* The bad fit in these 
variables is not important because their variations are small. Plots similar 
to those in figure 6 are presented in figure 7, where the parameters used were 
obtained by the ML method. The resulting residuals in the output variables are 
plotted in figure 8 . In this figure the standard errors of the output variables 
estimated from the residuals are also given. If the model of the airplane were 
correct, the residuals would form a white random sequence. Examination of their 
time histories in figure 8 indicates, however, the influence of certain modeling 
errors (incorrect form of aerodynamic model equations, uncorrected bias errors 
in measurement data, external disturbances). The incorrect model resulted in 
bias errors in the estimated stability and control parameters and in an increase 
of the Cramer-Rao lower bound. 
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The parameters obtained from all 16 runs mentioned are presented in fig- 
ures 9 and 10. Figure 9 includes the estimates by the EE method plotted against 
lift coefficient. The range of corresponds to the change in the angle of 

attack from approximately 0° to 11°. A distinction has been made between the 
resulting points where a pulse or doublet input was used. Each plot of a param- 
eter against C z was fitted by a first- or second-order polynomial, or by a com- 
bination of both. The change in the functional relationship from linear to qua- 
dratic was based on results from steady-state data as is shown later. 


In figure 10 the results from both estimation techniques were compared. 

The 2 a bounds on the EE estimates were computed from the differences between 
the estimates and the fitted curves. The 2a bounds can be used for an assess- 
ment of significant differences in results from both methods. As in table VII, 
the most pronounced difference exists in the estimates of c z a * The large dif- 


ferences in C 7 


and C Z( s; e at higher values of are not important because 


these parameters have small effect only on the resulting motion of the airplane. 
For the remaining parameters, both techniques provide, in general, equivalent 
estimates. Estimates of the parameters C Xf0 and are consistent even if 


the airspeed changes in the maneuvers analyzed were small. This is a promising 
indication of a possibility to estimate the performance characteristics of the 
airplane provided the necessary thrust information is available. 


Effect of nonlinear terms in aerodynamic model equations .- The examination 
of figures 9 and 10 shows variations of derivatives c z a ' c ma' and c n>q 

with C L at higher values of the lift coefficient. The effect of these varia- 
tions is not included in the aerodynamic model equations used. Therefore, for 
the assessment of model adequacy at higher angles of attack, the estimated stan- 
dard errors in the main output variables a v , q, and a z were plotted against 
Cl, as shown in figure 11. The estimates of measurement noise standard errors 
obtained from the compatibility check and instrumentation system characteristics, 
as presented in figure 3 and table II, have been added to these plots. The plots 
of figure 11 indicate the increase of the standard error in all three variables 
with increasing values of C^. A possible modeling error in the aerodynamic 
model equations has been checked using the data from one run at = 1.26. The 
data were analyzed by using the ML method and the modified aerodynamic model. 

The model included three additional nonlinear terms Cz a 2( a - ot Q ) 2 , 

qc (a - a 0 ) 

c m a 2( a _ a o) ' and Cmgq ~ ~* 


The resulting estimates are compared in table VIII with those based on the 
linear aerodynamics. In this table, the main diagonal elements of the informa- 
tion and parameter covariance matrices mjj and Sjj are also presented. The 
comparison shows that the extended model significantly influenced some of the 


parameters. The new estimates for C z and C 


Z S( 


are closer to those expected. 


With the nonlinear aerodynamics, there was an improvement in terms of standard 
errors in parameters, correlation between parameters, sensitivities, and measure- 
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ment noise standard errors. The values of nonlinear terms agree well with those 
predicted from figures 9 and 10. The standard errors of the nonlinear parameters 
were reasonably low, and any strong correlation between the linear and nonlinear 
parameters was not observed, as indicated by values of the correlation coeffi- 
cients given in table IX. 

Comparison of parameters estimated from steady-state and transient flight 
data . - The accuracy of estimated parameters was checked by the comparison of 
their values with those obtained from steady measurements. Data needed for the 
estimation of static stability and control derivatives were measured in slow 
acceleration-deceleration horizontal flights and in steady climbs and descents 
("saw-tooth" flights). In figure 12 the measured and fitted lift coefficient is 
plotted against the angle of attack. The same measured data from steady climbs 
with full power and steady descent with idle power are compared with the previous 
measurements in figure 13. 

The measured and fitted stick-fixed trim curves (S e against C L ) are plot- 
ted in figure 14 for two airplane center-of -gravity positions. The effect of 
power setting on these data is demonstrated in figure 15. All measured steady- 

state data indicate changes in derivatives C Za and with increasing 

and also with different power setting. 

The data from figure 12 and figure 14 were used for computing aerodynamic 
derivatives from expressions developed for an airplane with a conventional hori- 
zontal tail in reference 11 and modified for an all-movable tail. The relation- 
ship between directly estimated quantities from measured steady-state data and 
aerodynamic derivatives are summarized in appendix B. The results obtained are 
presented in table X where they are compared with the estimates from transient 
data. These estimates are the mean values from the range of C L within which 
the derivatives are assumed to have constant values. 

The comparison shows that the estimates from steady flights agree very 

well with those obtained from transient data by using the ML method. Only the 
values for the derivative differ significantly. However, the value 

for this derivative from steady-state data represents an approximation only. 

The results in table X can then be considered an indication of good confidence 
in the ML estimates for all parameters. 

In addition to the comparison in table X, two derivatives C Za and 

estimated from steady-state and transient data were plotted against C L in 
figure 16. All three estimates show the same trend with the increased value 
of C L . However, the differences in values of C Za remain. 

Effect of different input forms and power settings on estimated parameters . - 
In addition to the 16 runs analyzed, runs with different input forms were used 
to estimate the parameters. Data from runs with full and idle power were also 
used. All of these additional runs were measured at C L « 0.57. The examples 
of various stabilator-def lection time histories are given in figure 17. The 
input A consists of a series of pulses, the input B1 includes either sharp 
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pulses or double pulses, and the input B2 is a combination of sharp pulses with 
a slowly varying stabilator deflection used for better excitation of the air- 
speed changes. 


The effect of pulse and doublet on the parameter estimates has already been 
shown in figures 9 and 10. Some differences are apparent in derivatives c z a ' 


C^, and 



at higher angles of attack. 


The effect of the remaining input 


forms is shown in figure 18, where the new sets of estimates are compared with 
previous results. The results indicate that the estimates from runs with input A 
are very similar to those based on measurements with pulse and doublet input 
form. However, for the input forms B1 and B2, some differences in the estimates 
of pitching-moment derivatives appeared. These differences are mainly visible in 
the results of the EE method. For all input forms used, the Crame'r-Rao bound on 
the standard errors of the estimated parameters remained almost the same. Also, 
the strong correlation between C z and c zg e and between C m and c mg e was 


unchanged. These results agree with the investigation of effects of control 
inputs on the estimated parameters made in reference 7. 


For a more detailed investigation of a possible effect of the sharp pulse 
or short doublet on the estimates, two sets of results were compared. The input 
forms used and their harmonic contents are given in figure 19; the parameter 
estimates, relative sensitivities, standard errors, and the measurement noise 
standard errors are presented in table XI. 


The results from the run with sharp doublet show the incorrect sign in 



significant differences in C Za , and Cm ; decrease< 3 sensitivities; 

and increased inaccuracy. There is also a large standard error in the pitch- 
ing velocity. The time history of the residuals in q for this run included 
a dominant deterministic component whose form was similar to the expected time 
history of the angular acceleration in pitch. The degradation of results from 
the data with sharp doublet input can be attributed to modeling errors in the 
equations of motion and/or to the uncorrected errors in measured pitching veloc- 
ity. The comparison of measured and predicted frequency response curves dis- 
played in figure 20 shows this effect even more clearly. The predicted fre- 
quency response curves were computed from linearized equations (A2) , (A3) , and 
(A7) with the derivatives from table XI. The measured data were obtained by 
applying Fourier transforms to the measured time histories of 6 e , q, and a z . 


The measured data with slow doublet agree well with the prediction. In 
the second run with a sharp doublet, the frequency response curves a z ( joi) /6 e ( jw) 
again agree, but there are greater differences in the amplitude and phase charac- 
teristics of the frequency response function q( jto) /6 e ( jw) . The change in phase 
angle in the measured data is greater than theoretically possible considering 
the form of transfer function developed from equations (A2) and (A3) . To explain 
the differences mentioned, a new check and dynamic calibration of the pertinent 
rate gyro were made, and the measured pitch angle was compared with that pre- 
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dieted from the pitch rate gyro reading. Both checks indicated that the instru- 
mentation was functioning correctly and that no time-delay correction in mea- 
sured q(t) was needed. Another reason for the effect of the sharp doublet on 
the parameter estimates is in the modeling of the airplane. The confirmation of 
this conclusion, however, would need theoretical and experimental study which 
is beyond the scope of the present report. 

The results in figure 18 concerning power effect indicate that there were 
no differences between the estimates with full power and power required for the 
level flights. The changes in parameters with idle power agree with trends shown 
in measured steady-state data in figures 13 and 15, e.g., the decrease (in abso- 
lute value) of Cz a , Cm 5 e ' an ^ c mq an< ^ th e increase of C^. 


Lateral Characteristics 

For the estimates of lateral parameters, 28 maneuvers initiated from steady- 
state level flights at different airspeeds were available. The measured data 
were obtained from two flights. The results from eight repeated measurements at 
the same airspeed are given in table XII. They include the ensemble mean values 
and standard errors of the EE and ML estimates, and the average standard errors 
of a single measurement. The mean values obtained from results of both methods 
are quite consistent; the only significant differences were found in two less 
important derivatives Cyp and c l$ r ’ 

The ensemble standard errors of the EE estimates are, for most parameters, 
smaller than those of the ML estimates. In both cases the standard errors esti- 
mated from the ensemble do not agree with the standard errors of a single mea- 
surement. The ratio of these two different estimates varies between 1 and 5 for 
the EE method, and between 2 and 18 for the ML estimates. High values and vari- 
ability in these ratios could be caused by bias errors in the estimates resulting 
from various modeling errors and also from a small sample size. 

As an example of the comparison between measured and computed data, the 
time histories of one of the eight runs analyzed are presented in figures 21 
and 22. In figure 21, the computed responses are based on the EE estimates and 
in figure 22, on the ML estimates. These two figures indicate no significant 
differences for both sets of parameters, which is in agreement with results in 
table XII. 

The residuals from figure 22 and their standard errors are given in 
figure 23. In all time histories of residuals, some deterministic components 
are visible, the most significant of which are in the yawing velocity residuals. 
This study did not determine which modeling errors could contribute to the deter- 
ministic components in residuals. Because of the small amplitude of residuals 
in the given case, and also in the remaining runs, the existing modeling errors 
did not significantly influence the accuracy of the estimates. 

The estimated parameters from all 28 runs are plotted against C L in 
figure 24 and figure 25. In figure 24, the ML estimates were fitted by linear 
or quadratic polynomials. In figure 25, the ML estimates are compared with 
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those obtained by the EE method. As with the results in table XII, the dif- 
ferences in both estimates are not, in general, significant; furthermore, the 
results from the two flights do not differ substantially. 

The fitted curves in figure 24 can be used for the prediction of the lat- 
eral stability and control derivatives of the airplane. Table XIII shows the 
predicted values for C L = 0.62 and the standard error boundaries on the pre- 
diction errors. These boundaries represent the maximum and minimum standard 
error of the fitted curve © = 0(C L ) within the given range of C L . For com- 
parison, the last column of table XIII presents the Cramer-Rao lower bound on 
the estimated standard error of a single measurement. 

The estimates of the yawing-moment coefficient derivatives shown in fig- 
ures 24 and 25 revealed increased scatter with the increasing value of Cl- 
For the investigation of this trend, the standard errors of all output vari- 
ables were first plotted against Cl as shown in figure 26. However, these 
plots show only a moderate increase of this error for higher values of Cl- 
This increase indicates that the modeling errors do not change significantly 
for runs at higher angles of attack. 

The next step involved the investigation of sensitivities and standard 
errors in the parameter estimates. These characteristics deteriorate with 
increasing values of Cl which might lead to the conclusion that the input form 
used was not suitable for excitation at higher angles of attack. All the input 
forms used resulted in strong correlations between some of the stability and 
control derivatives, mainly between c 2p' c 2(5 e anc ^ ^p ' a % 7116 resu -*- ts i- n 

figure 26 demonstrate the increase in standard errors caused by modeling errors 
when one compares the estimates based on the equations of motion, on kinematic 
equations only (compatibility check) , and on the instrumentation system alone. 

The difference in the standard errors s(r) in two flights was probably caused 
by more pronounced modeling errors in flight 21 , where the input was similar to 
that in flight 26, but the rudder amplitude was higher. 

Comparison of results from transient and steady-state measurements .- As 
with the longitudinal case, the results from the lateral transient flight data 
were compared with those obtained from steady nonsymmetric flights (steady 
straight sideslips) . Unfortunately, the lateral steady-state data cannot pro- 
vide estimates of aerodynamic derivatives directly without additional a priori 
information. For that reason, the measured relationships between the sideslip 
angle and the bank angle, the aileron deflection and the rudder deflection were 
compared with those predicted using the parameters estimated from transient 
data. If linearized lateral equations of motion are assumed, it can be shown 
that the aforementioned relationships are linear and that their slopes and inter- 
sects depend on certain combinations of derivatives. The analytical form of the 
relationships 6 = 8(4>), 8 = 8 (S a ) and 8 = 8(6 r ) is presented in appendix B. 

In figures 27 to 29, the measured lateral static aerodynamic characteristics 
for two values of C L are presented. The estimates of those characteristics 
were computed from equations (B13) to (B15) using the ML estimates for the param- 
eters. The agreement between measurements and computed lines is generally good. 
The only discrepancy is for the intersection of the line 8(<5 r ) at Cl = 1.19. 
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This intersection depends primarily on the parameter C nf0 . The estimates of 
this parameter were incorrect as to sign. One explanation could be the exis- 
tence of modeling errors in the yawing-moment equation; however, the proof would 
require more detailed analysis. 

The results in figures 27 to 29 verified to some extent the estimates of 
the derivative c yg » and the combination of derivatives c 2 < 5 a / c 2|3 and C n<5 r / C ng* 

For the verification of the control derivatives c 2g a and C n g r , the aileron and 

rudder effectiveness were measured in steady flights with 6 ~ 0 and with the 
additional rolling and yawing moments provided by the onboard rocket system. The 
resulting derivatives from these measurements are given in figure 30 and compared 
with the ML estimates taken from figure 24. The steady-state data resulted in 
slightly higher absolute values for both control derivatives than for the tran- 
sient data. The differences are, however, not significant. 

Effect of different input forms and power settings on estimated parameter s .- 
All previous lateral data analyzed included only one type of input form for both 
the rudder and aileron deflection. To investigate for an effect of other input 
forms, five other types of input forms shown in figure 31 have been used. The 
input A1 is similar to the standard input, but the aileron doublet is shorter and 
has greater amplitude. The input A2 is composed of input A1 repeated three times 
In the input B1 the rudder pulse is followed by the aileron doublet. As in the 
previous case the input B2 is composed of input Bl repeated three times. Finally 
in the input C the aileron doublet precedes the rudder doublet. 

Estimated parameters from several runs with the inputs described are pre- 
sented in figure 32. They are compared with the ensemble mean values given in 
table XII which were estimated from runs where the standard inputs were used. 

In all cases the parameters are the ML estimates. 

The results in figure 32 revealed a significant effect of input forms on 
virtually all estimated parameters. The changes in the estimates are, in gen- 
eral, greater in the derivatives of Cy and C\ than those in derivatives of 
C n if the ensemble mean value and its confidence interval are considered as a 
reference. In many cases it was not possible to determine which estimates are 
more accurate because the standard errors and the correlation coefficients were 
almost the same for significantly different values of the same parameter. 

The variability in the lateral parameter estimates due to different input 
forms was found to be less pronounced in the study covered by reference 7. This 
could be explained by the dependence of the input form sensitivity on the char- 
acteristics of an airplane. The impossibility of selecting the best estimates 
of parameters from runs with different inputs demonstrates the existing problem 
of an optimal input form for parameter identification and the problem of accu- 
racy assessment of the parameter estimates. 

In figure 32, the effect of different power settings on the parameter esti- 
mates is also shown. The results from runs excited from steady-state climbs 
with full power agree with those from the measurements in horizontal flight. 

This agreement was expected because in both experiments the power setting dif- 
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fered only slightly. On the other hand, the results from runs with idle power 
show a decrease in the absolute value of static derivatives Cy^ and c ng' and 

control derivatives c Yg r and c n< 5 r - These changes and their directions have 

been expected. They are caused by the slipstream effect on the sidewash angle 
and the dynamic pressure at the tail. 


ASSESSMENT OF ESTIMATED PARAMETERS 

A comparison of the results in table X shows the ML parameter estimates 
to be the best values for the longitudinal stability and control derivatives 
of the airplane under test within the range of C L where these derivatives have 
constant values, e.g., < 0.9. Power will not significantly change the value 

of these derivatives, with the exception of power settings close to zero-power 
conditions. For both methods, the standard error of all important derivatives 
is about 2 percent. This error can also be a measure of the overall accuracy 
because of the agreement between the results from steady-state and transient 
data, and between the EE and ML estimates as shown in figure 18 for the pulse 
input doublet and combination of both (input A) . The values of derivatives for 
C L > 0.9 can be obtained from figure 10. The accuracy of these derivatives 
will deteriorate with the increasing value of because the estimates are 

influenced by the input form, power conditions, and the uncertainty in the aero- 
dynamic model equations. The best values of the stability and control deriva- 
tives for the lateral motion are given by the fitted curves in figure 24. The 
bounds on their standard errors can be obtained from table XIII. The standard 
errors for the aerodynamic static derivatives are between 1 and 3 percent, for 
the damping derivatives between 2 and 10 percent, and for the primary control 
derivatives between 2 and 5 percent. For some less significant derivatives the 
standard error can be as high as 70 percent. There is a good agreement between 
some derivatives and their combinations obtained from steady-state and transient 
data. 


The EE and ML estimates of all important derivates also agree provided 
that the data with the same input form were used in the analysis. On the other 
hand, the estimated derivatives depend strongly on the input form used which 
means that the input-form dependence degrades the accuracy of the estimates. 

The effect of power is insignificant with the exception of the idle power 
regimes . 


CONCLUDING REMARKS 

A complete set of stability and control parameters included in the aero- 
dynamic model equations was estimated from maneuvering flight data for a low- 
wing, single-engine, general aviation airplane. Most of the estimated param- 
eters obtained by using the equation error and maximum likelihood methods 
agreed within two-standard-deviation confidence intervals for the parameter. 
This agreement was made possible because sufficient accuracy of measured data 
was achieved from a thorough ground and flight calibration of the instrumenta- 
tion system, a check of the system before and after each flight, and correction 
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of the measured data for bias errors determined from compatibility checks of 
measured response variables. The estimated static parameters also agreed with 
the results from steady flights. The comparison of parameters obtained from 
different methods and through repeated measurements resulted in the determina- 
tion of the best values for the estimated parameter and in the specification 
of their accuracies. 


Longitudinal Parameters 

Both estimation methods provided identical values for most of the param- 
eters. The significant difference was found in the derivative of the vertical- 
force coefficient with respect to the angle of attack. The maximum likelihood 
estimates agreed better with the computed parameters from steady-flight data than 
with those from the equation error method. The engine power did not change the 
values of parameters significantly, with the exception of power settings close 
to zero-power conditions. The standard error of all important parameters was 
about 2 percent. The accuracy of parameters deteriorated with the increasing 
value of the lift coefficient. The high angle of attack and rapid maneuvers 
created some uncertainties in aerodynamic model equations. The addition of non- 
linear aerodynamic terms could improve the parameter estimates at high values 
of the lift coefficient. 


Lateral Parameters 

The results from eight repeated measurements under the same flight condi- 
tions showed that the mean values from both estimation methods agreed in general. 
The only significant differences were found in two less important derivatives. 

The ensemble standard errors of parameters obtained by the use of the equation 
error method were smaller than those of the maximum likelihood estimates. At 
the same time, the ensemble standard errors of the parameters from both methods 
were higher than the standard errors of a single measurement. 

The comparison of parameters plotted against the lift coefficient indicated 
no significant differences between the results of the two methods. The standard 
errors of the main derivatives varied between 1 and 10 percent. For some less 
important derivatives, this error was as high as 70 percent. There was good 
agreement between combinations of the static stability and control parameters 
obtained from steady and transient data. The estimated parameters depended 
strongly on the input form used. The effect of power was not significant with 
the exception of the idle power regimes. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 26, 1979 
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APPENDIX A 


AIRPLANE STATE AND OUTPUT EQUATIONS 

The airplane equations of motion are referred to the body axes. They are 
based on the following assumptions: 

(1) The airplane is a rigid body. 

(2) The stabilator deflection excites only the longitudinal motion whereas 

the rudder and ailerons excite the lateral motion which is coupled 
with the longitudinal motion. 

(3) There are no external disturbances to the airplane. 

(4) The aerodynamic model equations include also the effect of propeller 

and have the form 


C X = c X,o + c X a 

C Y = c Y,o + c Yg(3 ~ 6 0 ) + C Yp ^ + C Ygr (S r - 6 r ,c>) 

qc 

C Z “ c Z,o + “ «o) + c Z q ~ + c Z 5 e ( 6 e “ 5 e,o> 


c l = c l,o + C 2 B (B - B 0 ) + c z p ^ + C lr — + C Zga (5 a - 6 a ,o> 
+ c l6 r < 5 r - Sr ,o) 

ac qc 

Cm = ^,0 + c m a ( a “ a o) + C^ ~ + c mq ~ + '-mgg^e “ ^e,o^ 

c n = c n,o + c n B (B - 3 0 ) + c n p ~ + c n r ~ + c n 5 a ^a ~ ^a,o) 
+ ^ n 6r ~ < "* r /o) 
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APPENDIX A 


where 


C x = C»ji cos a T + Cl sin a - Cq cos ol 


C z = -Cip sin a T - C L cos a - c d sin a 

and where the index o denotes the value of the coefficient, output, or input 
variable corresponding to the initial steady-state flight conditions. 

These assumptions allow the longitudinal state equations to be expressed 
as 


. qS r— 

u = -qw - g sin 0 + — |C X , 0 + C Xa (a - a Q )J 


(AT) 


. qS 

w = qu + g cos 0 + — 

m 


qc 


'Z ,o + c Z a ( a - a o> + c Z q ~ + c Z 5 e ( 5 e " 5 e,o) 


(A2) 


. qSc 

q 'i7 


qc 

c m,o + c i% ( a _ a o) + ~ + c m§ e ^e ~ ^e,o) 


(A3) 


0 = 


(A4) 


and the output equations to be expressed as 


a v = tan 


-1 


w - qxct' 


1 . 

a x =-(u+qw+g sin 0) 

g 


1 . 

az = ~(w - qu - g cos 0) 
9 


(A5) 


(A6) 


(A7) 
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In the pitching-moment equation 


^m / o “ ‘-iTi/O + c m, 


( pSc gc 

Cz '° + ^ 


cos e 0 


pSc 

' m a “ " m a ^ 17 Lm 7 z ct 


°Aq 


- % + 



pSc 

Sse " Sse + 7T Cm a Cz 6e 


True airspeed and angle of attack are computed from the equations 


V = 



+ w 


2 


a 


tan 



For the equation error method the state equations were modified as 


gm 

=1 a X = c X,o + c X a (a - a 0 ) 


gm qc 

^ a z = c z,o + c z a (a - a Q ) + C Zq — + c z^ e (5 e - 6 e ,o) 


fi q 2v <$ e 


l y . 


qc 


5^ q - + C "b (a ao> + C ™q 2V + Sse^e " 6 e,o> 


(A8) 


(A9) 


(AT 0) 


(All) 


(AT 2 ) 


(A13) 


(AT 4) 
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The lateral equations can be expressed as 


v 


= -u E r + w E p + g cos 


0 E 


sin 4> + 



+ C Y g (B - 0 O > 


+ C\ 


pb 


L P 2V 


w + C Y6r (6 i 



P = 


ft 


I X 1 Z - I xz‘ i 


+ : F 2 

I X I Z “ *xz 


r = 


I X I Z “ I XZ 2 


T XZ 

F2 + ; ft 

¥z ~ I xz 


<£> = p + (q E sin 4> + r cos 4>) tan 0 E 


where 


F-| = (I Y ~ ~ i XZ < JeP + < i sb 


C 1 ,o + C Iq($ ~ ^o) + c l p 2V 


rb 


+ c l r — + c 2s a ( 6 a - 6 a,o> + c 25 r (6 r 5 r,o) 


F2 = (I x ~ ly^EP " I XZ < 3E r + 3 sb 


c n,o + c ne ( ® ^o) + c n p 2v 


+ c n r ~ + c n<$ a ( 6 a “ 5 a,o> + c n< 5 r ( 6 r 6 r,o) 


V = 



+ w 


2 

E 


3 = sin 1 


v 

V 


(A15) 


(A16) 


(A17) 


(AT 8) 


(AT 9) 


(A20) 
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u E « V E cos a VE cos 8 Ve 

W E « V E sin a VE cos 8 Ve + q E *a - 

the index E denotes the measured quantity. 

The output equations related to the lateral motion have the form 


8 V = sin 


-1 


v + r 


xg - pzg ' 


V 


1 . 

a y = -(v + u E r - wsp - g cos 0 E sin <j>) 

g 


For the equation error method, the lateral equations were modified 


gm pb 

— a Y = c Yf0 + c Y g (8 - 8 0 ) + c Yp — + c Yfir (6 r - 6 rf0 ) 


ix 

qSb 



(l Y ~ Iz\ 


I XZ 

(pq + r) 

ix 


- C2 ,o + c ?8 - e o) + % % + C ?r Tv + C Z5a (6 a " 6 a,o> 

+ C 2(Sr (6 r - 5 r , 0 ) 

I Z 

qSb 

= c n,o + c ng - 8 0 ) + c n p “ + c n r ~ + Cngg^a " 5 a,o) 
+ Cn 6r ^ r - ^ r '° ^ 



(A21) 


(A22) 


as 


(A23) 


(A24) 


(A25) 
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RELATIONSHIP BETWEEN AERODYNAMIC DERIVATIVES AND PARAMETERS 
ESTIMATED FROM STEADY-STATE DATA 
From the measured steady-state data of the aircraft under test 

C L = CL(ot/h,6t»P ower setting) 

6 e = S e (C L ,h,<$ t/ power setting) 

the following parameters for given power and trim tab setting can be determined: 

Slope of the Cj^a) curve, dC^/da 
Slope of the <5 s (Cl) curve, dfig/dC^ 


VrpS i 


Ah 



e 


(Bl) 


and 


Ah dS e 



(B2) 


where 


1 s t z 

Vrp (B3) 

1 + F Sc 


ch no 

F = (B4) 

*n 
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a l 


3C 


L,t 


3a t 


and where AS e is an increment due to change of relative center-of -gravity 
position Ah for a given C L . 

Using these parameters allows the static stability and control derivatives 
to be obtained from the expressions 


C ®6« 


-v T a l 


1 + — H r 


(B5) 


CL Se “ i SSe 
1 1 


(B6) 


dC L /da 


dS e 

1 + 35 


(B7) 


c % = ~ Cl rx Hn 


(B8) 


The tail contribution to the damping derivatives can be approximated as 


S t /] t \2 l t 

= -2a-| —I — I - -2V Tai — 
*3/ 1 S \c / c 


(B9) 


1 1 3e 

^ t " _2VTai r ^ 


(BIO) 


c Z q )t “ -2^1 


(Bll) 


where the rate of downwash angle can be approximated as 


28 


APPENDIX B 


9e 2Cl cx 

1 ' at — — 

9a TTA 


(B1 2) 


For comparison of the estimated aerodyanmic derivatives from steady-state 
and transient flight data it was assumed that 


Cz a ** CL a 
Cz 6e “ " C L6e 



The modified derivatives and Cj{^ were computed from equations (A9) 

and (A10) . In the given case for he = 0.206 it was found that Cj^ = -1 .015, 
whereas Cj^ = -0.80. However, the difference between c m 5 e and c mg e was 
negligible. 

In steady-state nonsymmetric flights the following relationships can be 
measured : 


3 = 3 ( 4 >) 

6 = 3 ( 6 a ) 
6 = B(S r ) 


The analytical ^ forms for these relationships are developed from equations (A15) 
to ( A 1 7) for p = r= p= q = r= 0 and B 0 = 6 af0 = <$ r fQ = 0 as 


K <t> / Y fi r c n,o k 5a c 2 ,o\ 

C Yg \ ' c n 5 r 1 “ k 6a k 6r / c Yg 


(B13) 


l,o k Sr^n,o ^fia 1 ~ k <Sa k 6r 

3 <5 a (B14) 

c 2g(1 - k 6r k n2B^ c 2g 1 ~ k <5r k n2B 


^n ,o k 6a^2 ,o ^ n <5r ^ k <5a k Sr 

B = - 6 r (B1 5) 

C n g(1 k 6a k 2nB^ c ng 1 ~ k 6a k 2 nB 
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In these equations. 


2mg 

C 7 = - cos 0 COS <t> 

psv 2 


= 


1 - 


C *Sr C ng 1 - k 6a k Zne 
^n( 5 r ^ ~ k 5a k 6r 


k <5a 


n 6< 

C *<5< 


k Sr 


Cl 6r 

Cn 6r 


k Zn£ 


1 

k nZ6 


c Zg 

c ng 


(B16) 


(B 1 7) 


(B18) 


(B19) 


Equations (B13) to (B15) also define the relationship between the slope 
and intersect of the lateral steady-state characteristics and the lateral sta- 
bility and control derivatives. 
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TABLE I.- GEOMETRIC, MASS, AND INERTIA CHARACTERISTICS OF AIRPLANE 


Wing: 

Area, m 2 13.56 

Aspect ratio 7.35 

Span, m 9.98 

Mean aerodynamic chord, m 1.34 

Ailerons : 

Area, m 2 0.94 

Half-span, m 1.65 

Vertical tail: 

Area, m 2 1.36 

Rudder area, m 2 0.43 

Horizontal tail (stabilator) : 

Area, m 2 2.51 

Aspect ratio 4.21 

Span, m 3.25 

Tail length (c.g. position at 0.206c), m 4.21 

Fuselage length, m 7.85 


Mass: 

Aircraft mass at take-off, kg 974 

Aircraft mass at landing (no fuel) , kg 877 

No fuel Full fuel 

Inertia: 

l x , kg-m 2 1568 1888 

l y , kg-m 2 2125 2142 

l z , kg-m 2 2326 3557 

I xz , kg-m 2 140 142 
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TABLE II.- CHARACTERISTICS OF THE INSTRUMENTATION SYSTEM 


Quantity measured 

Transducer 

Range 

(a) 

Static 

sensitivity 

(b) 

Resolution 

(c) 

rms measurement 
error 
(c) 

Unit 

(b) 

Percent of 
full range 

Longitudinal acceleration, g units 


-1 

to 

1 

2.54 

0.001 

0.0046 

0.23 

Lateral acceleration, g units 

l Servo accelerometer 

-1 

to 

1 

2.48 

.001 

.0050 

.25 

Vertical acceleration, g units 

J 

-3 

to 

6 

.56 

.001 

.0050 

.06 

Rolling velocity, deg/sec 

> 

-102 

to 

102 

.025 

.12 

.20 

.10 

Pitching velocity, deg/sec 

l Rate gyro 

-29 

to 

29 

.088 

.032 

.19 

.33 

Yawing velocity, deg/sec 

J 

-29 

to 

29 

.084 

.034 

.080 

.14 

Roll angle, deg 

'S Vertical gyro 

-90 

to 

90 

.028 

.10 

.077 

.04 

Pitch angle, deg 

J 

-87 

to 

87 

.029 

.098 

.092 

.05 

Angle of sideslip, deg 

~\ Flow direction 

-12 

to 

27 

.127 

.029 

.027 

.07 

Angle of attack, deg 

velocity sensor 

-29 

to 

32 

.124 

.018 

.019 

.03 

Right aileron angle, deg 


-23 

to 

10 

.147 

.0020 

.019 

.06 

Left aileron angle, deg 

l Control position 

-10 

to 

25 

.142 

.0020 

.0061 

.02 

Stabilator angle, deg 

/ transducer 

-16 

to 

3 

.263 

.010 

.0037 

.02 

Rudder angle, deg 

J 

-31 

to 

28 

.084 

.011 

.0091 

.02 

Airspeed, m/sec 

Pressure transducer 

0 

to 

75 

.067 

.037 

.89 

1 .2 

Altitude, m 

Altimeter 

-150 

to 

2900 

.0016 



— 

Air temperature, °C 

Thermometer 

. 

-18 

to 

38 

.089 

- 


— 


a Working range of the channel. 

Obtained as volts per pertinent unit. 

c Ref erred to a reading from the digitized tape 



TABLE III.- DYNAMIC CHARACTERISTICS OF INSTRUMENTATION SYSTEM 


Quantity measured 

Natural 

frequency, 

Hz 

Damping 

ratio 

Equivalent 

time 

constant, 

sec 

Longitudinal acceleration, g units 

402 

1 .58 

0.0012 

Lateral acceleration, g units 

216 

1.10 

.0016 

Vertical acceleration, g units 

921 

1.58 

.0005 

Rolling velocity, deg/sec 
Pitching velocity, deg/sec 

1 27 

.64 

.0075 

Yawing velocity, deg/sec 
Angle of sideslip, deg 

J 

(a) 

(a) 


Angle of attack, deg 

23 

.085 

.0012 

Airspeed, m/sec 

— 




a At V = 50 m/sec. 
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TABLE IV.- FLIGHT CONDITIONS AND AVERAGE MASS AND INERTIA CHARACTERISTICS 

OF AIRPLANE IN TEST FLIGHTS 


Flight m ' h ' _ I X» 

kg percent c kg-nr 


r XZ' 

kg-m z 

P# 

kg/m^ 

Experiment 

140 

1 .076 

Longitudinal responses to short duration 
pulses; lateral responses; v^ ~ 85 knots 

130 

1 .107 

Lateral responses; 62 knots < V^ < 100 knots 

— 

1 .044 

Longitudinal responses; 

62 knots < V^ < 100 knots 

130 

— j 

• 

o 

CO 

OJ 

Lateral responses; V^ * 85 knots; different 
input forms and power settings 

— 

1 .066 

Longitudinal responses; Vj •« 85 knots; 
different power settings 

— 

.989 









TABLE V.- ESTIMATED INSTRUMENT BIAS ERRORS FROM STEADY- STATE AND 
TRANSIENT DATA FOR LONGITUDINAL MOTION 


Parameter 

Mean value 
from steady- 

Estimate from 
(« 

transient data 
i) 


state data 

No fixed value 

One fixed value 

b a X ' g units 

0 

0.000 

(0.0030) 

b 0 

b a z » 9 units 

-.035 

-.022 

( .0067) 

-.025 (0.0067) 

bq, deg/sec 

.068 

-.01 

( .042 ) 

-.01 ( .042 ) 

bar deg 

-.14 

-.28 

( .078 ) 

-.28 ( .080 ) 

bg, deg 

.81 

.77 

( .063 ) 

.74 ( .063 ) 

X a 

0 

.01 

( .0030) 

.01 ( .0030) 


a Numbers in parentheses are Cramer -Rao lower bounds on standard 
errors. 

b Fixed value. 









TABLE VI.- ESTIMATED INSTRUMENT BIAS ERRORS FROM STEADY-STATE AND 
TRANSIENT DATA FOR LATERAL MOTION 


Parameter 

Mean value 
from steady- 
state data 

Estimate from transient data 
(a) 

No fixed value 

One fixed value 

b a Y ' g units 

0 

-0.017 (0.0040) 

b 0 

bp, deg/sec 

3.58 

2.69 ( .042 ) 

2.69 (0.042 ) 

b r , deg/sec 

.55 

.24 ( .037 ) 

.26 ( .036 ) 

bg, deg 

0 

.00 ( .057 ) 

-.02 ( .057 ) 

btyf <3eg 

1 .59 

.88 ( .080 ) 

.09 ( .080 ) 

*8 

0 

.008 ( .0021) 

.003 ( .0021) 


a Numbers in parentheses are Cramer -Rao lower bounds on standard 
errors . 

b Fixed value. 
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TABLE VIII.- PARAMETERS ESTIMATED PROM TRANSIENT FLIGHT DATA USING TWO 


DIFFERENT AERODYNAMIC MODELS AND THE MAXIMUM LIKELIHOOD METHOD 
[Flight 25, run 79B; C L = 1.26] 


Parameter 

Linear aerodynamics 

Nonlinear aerodynamics 

Estimate, 

© 

Sensitivity, 

1/mjj, 

percent 

Standard 
error , 
Sjj, 
percent 

Estimate , 

A 

© 

Sensitivity, 
1 / m j j r 
percent 

Standard 
error , 
Sjj, 
percent 

0 

X 

a 

0.252 

0.17 

0.31 

0.251 

0.11 

0.25 

Cx a 

1.15 

1 .0 

1 .5 

1 .17 

1 .1 

1 .5 

c z,o 

-1 .214 

.09 

.25 

-1 .238 

.04 

.20 

c z a 

-3.44 

1 .3 

2.1 

-3.22 

.80 

2.2 

c z a 2 


— 

— 

1 1 

2.0 

9.6 

% 

-33 

3.1 

6.8 

-23 

3.6 

7.4 

C Z6e 

-2.67 

3.9 

8.7 

-2.08 

3.0 

8.2 

*-^1,0 

-.0097 

3.1 

6.6 


3.5 

17 


-1 .301 

.65 

1 .0 

-1 .473 

.21 

1 .1 

C i% 2 


— 

— 

-4.8 

.66 

5.4 


-18.6 

.97 

2.4 

-15.6 

.40 

2.5 



— 

— 

52 

2.6 

14 

°ili6e 

-3.44 

.56 

1.4 

-3.12 

.18 

.97 

s (V) , m/sec 

.19 

— 

— 

.17 

— 


s (a) , deg 

.33 

— 

— 

.24 

— 


s(q) , deg/sec 

.77 

— 

— 

.61 

— 


s(0) , deg 

.43 

— 

— 

.12 

— 


s(ax) , 9 units 

.0049 

— 

— 

.0057 

— 


s(a z ) , g units 

.087 

— 

— 

.029 

— 
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TABLE X.- LONGITUDINAL PARAMETERS ESTIMATED FROM STEADY-STATE 
AND TRANSIENT FLIGHT DATA 


Der ivative 

Steady-state 

data 

Transient data 
___ (a) 

EE method 

ML method 

b c z 
z a 

b -5 . 1 0 

b -4 . 68 

(0.04) 

b -5 . 3 

(0.1 ) 

c z 

z q 

C -1 8 

-16 

0 ) 

b_i9 

(3 ) 

Cz 6e 

-1.04 

-1.02 

( .06) 

-1.2 

( .2 ) 

d r> * 
Lm a 

d -.80 

d - .68 

( .03) 

d - . 80 

( .02) 

d c 1 

'-Trig 

c -30 . 6 

d -27 . 3 

( .2 ) 

d -24 . 2 

( -4 ) 

Cmfie 

-3.26 

-3.32 

( .04) 

-3.32 

( .05) 


a Numbers in parentheses are standard errors of ensemble mean. 

b For C L ^ 1.0. 

c Tail contribution only. 

d For C L S 0.7. 
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TABLE XI.- PARAMETERS ESTIMATED FROM TRANSIENT FLIGHT DATA WITH 


TWO DIFFERENT INPUTS BY MAXIMUM LIKELIHOOD METHOD 



a Estimates for he = 0.206. 


TABLE XII.- PARAMETERS AND THEIR STANDARD ERRORS ESTIMATED FROM REPEATED MEASUREMENTS 


USING TWO ESTIMATION METHODS 



a Ensemble mean value. 


^Ensemble standard error . 


c Average standard error of estimates. 

CaJ 

















.60 | — 


s(V) = 0.089 m/sec 





Figure 3.- Time histories and standard errors of residuals. 
Longitudinal motion; flight 25, run 13B. 
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Figure 4.- Time histories of measured and predicted output variables 

Lateral motion; flight 21, run 26. 
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Figure 7.- Measur 
computed by us 
Flight 25, run 




















EE method: 


ML method: 


Mean value O Pulse input 










Figure 12.- Measured and fitted lift coefficient plotted against angle of 
attack. Acceleration-deceleration levels. 
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Figure 16.- Comparison of two parameters estimated from steady-state and 

transient flight data. 












Mean value. Pulse and doublet 
input, level flights 
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Open symbols: EE estimates 

Closed symbols: ML estimates 
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Figure 18.- Effect of different input forms and power settings on estimated 

longitudinal parameters. 
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□ Mean value. Pulse and doublet 
input, level flights 

J 2 a bound for estimates 



Open symbols: EE estimates 

Closed symbols: ML estimates 



A B1 B2 Full Idle 

Input Power 


Figure 18.- Concluded 
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deg 



Figure 20.- Comparison of measured frequency response curves for two different 
inputs with those computed by using parameters obtained by maximum likeli- 
hood method in time domain. 
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ted 




and those computed by 
Flight 21, run 26. 
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Figure 22.- Measured lateral 
using parameters obtained 


flight data time histories and 
by maximum likelihood method. 


those computed by 
Flight 21 , run 26. 
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o Flight 21 

• Flight 26 



0 .2 .4 .6 .8 1.0 1.2 


Figure 24.- Estimated lateral parameters from flight data. Maximum 

likelihood method. 
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Figure 25.- Comparison of lateral parameter estimates obtained from flight 
data using equation error and maximum likelihood methods. 
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Figure 31 Different aileron and rudder deflections used in test 




□ Mean values 
I 2a bound for estimates 
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A1 A2 B1 B2 C Full Idle 

Input form Power 

Figure 32.- Effect of different input forms and power settings on estimated 
lateral parameters. Maximum likelihood method. 
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□ Mean values 
X 2 a bound for estimates 





A2 B1 B2 
Input form 


Figure 32.- Continued. 
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□ Mean values 
I 2a bound for estimates 



A1 A2 B1 B2 C Full Idle 

Input form Power 


Figure 32.- Continued. 
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□ Mean values 
X 2a bound for estimates 


Input form 


Full Idle 
Power 


Figure 32.- Concluded. 
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